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Abstract: We calculate anomaly induced conductivities from a holographic gauge 
theory model using Kubo formulas, making a clear conceptual distinction between 
O ! thermodynamic state variables such as chemical potentials and external background 

fields. This allows us to pinpoint ambiguities in previous holographic calculations 
of the chiral magnetic conductivity. We also calculate the corresponding anoma- 
lous current three-point functions in special kinematic regimes. We compare the 
holographic results to weak coupling calculations using both dimensional regulariza- 
tion and cutoff regularization. In order to reproduce the weak coupling results it is 
necessary to allow for singular holographic gauge field configurations when a chiral 
chemical potential is introduced for a chiral charge defined through a gauge invari- 
ant but non-conserved chiral density. We argue that this is appropriate for actually 
addressing charge separation due to the chiral magnetic effect. 
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1. Introduction 

The chiral anomaly of QED is responsible for two particularly interesting effects of 
strong magnetic fields in dense, strongly interacting matter as found in neutron stars 
or heavy-ion collisions. At large quark chemical potential fi, chirally restored quark 
matter gives rise to an axial current parallel to the magnetic field [1-3] 

J 5 = ^B, (1.1) 

which may lead to observable effects in strongly magnetized neutron stars [4]. 

In heavy-ion collisions, one expects initial magnetic fields that momentarily ex- 
ceed even those found in magnetars. It has been proposed by Kharzeev et al. [5-9] 
that the analogous effect [10] 

J = ^ 5 B, (1.2) 

where J is the electromagnetic current and ^ a chemical potential for an asymmetry 
in the number of right and left chiral quarks, could render observable event-by-event 
P and CP violations from topologically nontrivial gluon configurations. Indeed, there 
is recent experimental evidence for this "chiral magnetic effect" (CME) in the form of 



- 1 - 



charge separation in heavy ion collisions with respect to the reaction plane [11,12] (see 
however [13, 14]), whose normal vector is expected to coincide with the direction of 
strong initial magnetic fields. For lattice studies of the effect, see for example [15,16], 

The anomalous conductivities (|1 . 1| ) and ( |1.2|) have recently also been studied 
in holographic models of QCD by introducing chemical potentials for left and right 
chiral quarks as boundary values for corresponding bulk gauge fields [17, 18]. How- 
ever, it was pointed out by Ref. [19] that in these calculations the axial anomaly was 
not realized in covariant form and that the corresponding electromagnetic current 
was not strictly conserved. Correcting the situation by means of Bardeen's counter- 
term [20,21] instead led to a vanishing 1 result for the electromagnetic current in the 
holographic QCD model due to Sakai and Sugimoto [23, 24] 2 , while recovering the 
result ( |1 . 1| ) for the anomalous axial conductivity. 



Indeed, the two anomalous conductivities ( |1.1|) and ( |1.2|) differ in that in the 
former case there is no difficulty with introducing a chemical potential for quark 
number, while a chemical potential for chirality refers to a chiral current that is either 
gauge invariant and anomalous or conserved but not gauge invariant. In Ref. [10], 
the chiral magnetic effect ( |1.2| ) was shown to be an exact result when the chiral 
chemical potential is conjugate to conserved chiral charges that are, however, only 
gauge invariant when integrated over all of space in spatially homogenous situations 
(this point was most recently also made in Ref. [26]). 

However, charge separation in heavy-ion collisions clearly calls for inhomogeneous 
situations, since with V ■ B = we have 

i = -V. J = -^B.V*. (!,) 

It would therefore appear important to introduce a chiral chemical potential con- 
jugate to gauge invariant axial currents, despite them being anomalous. At least 
as long as the electric field is zero and as long as the chiral charge decay rate is 
suppressed (as it is in the large- N c limit [27], and the fact that we indeed find time- 
independent solutions in the presence of a chiral chemical potential is a reflection 
of this fact), it should be admissible to consider chemical potentials defined with 
respect to the gauge invariant chiral density in a thermodynamic description. Such a 
chiral chemical potential thus serves as a model parameter for the imbalance between 
the number of left-handed and right-handed fermions that is assumed to be induced 
by topologically non-trivial gluon field configurations during the out-of-equilibrium 
early stages of a heavy ion collision. 



1 It has been argued that this is an artifact of the grand canonical ensemble and that the weak 
coupling result would be recovered in the canonical ensemble [22]. 

2 In Ref. [25] a Unite result was obtained in a bottom-up model that is nonzero only due to extra 
scalar fields. 
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It is however important to distinguish between thermodynamic state variables 
such as chemical potentials and background gauge fields (as also pointed out by Ref. 
[26]). The holographic dictionary instructs us to construct a functional of boundary 
fields and n-point functions are obtained by functional differentiation with respect 
to the boundary fields. For a gauge field, the expansion close to the boundary takes 
the form 

A,(x,r)=Af( X ) + ^^ + .... 

The leading term in this expansion is the source for the current J M . The sub- leading 
term is often identified with the one-point function of the current. This is, however, 
not true in general. As has been pointed out in Ref. [19], in the presence of a bulk 
Chern-Simons term, the current receives also contributions from the Chern-Simons 
term and Afr (x) can, in general, not be identified with the vev of the current. On 
the other hand, a constant value of A Q is often identified with a chemical potential. 
This is, however, slightly misleading since the holographic realization of the chemical 
potential is given by the potential difference between the boundary and the horizon 
[28] and only in a gauge in which A vanishes at the horizon such an identification 
can be made. Even in this case, we have to keep in mind that the boundary value of 
the gauge field is the source of the current whereas the potential difference between 
horizon and boundary is the chemical potential. We will keep this distinction explicit 
in this paper. 

In section 2 we shall show that by distinguishing the chemical potential from 
the background gauge fields one can reproduce the usual result ( |1.2j ) for the chiral 
magnetic effect when /j,$ refers to the gauge invariant chiral current 3 and for a strictly 
conserved electromagnetic current J . However, this requires singular gauge field 
configurations in the bulk of AdS space, which appears to be the price to pay for 
having introduced a chemical potential for an anomalous charge. We also reproduce 
the uncontested result ( |1 . 1|) for the axial current at finite quark chemical potential 
and magnetic fields, as well as a new anomalous conductivity, albeit one of perhaps 
mere academic interest as it refers to nonzero axial magnetic fields. Moreover, we 
derive results for anomalous three-point functions in certain kinematic limits. In 
section 3 we reproduce all these results in weak coupling calculations using gauge 
invariant dimensional regularization without and cutoff regularization with the need 
for introducing Bardeen's counter-term. 



2. Holographic Kubo formulas for anomalous conductivities 

We will consider the simplest possible holographic model for one quark flavor in a 



3 Thc chiral current is gauge invariant under the non-anomalous vectorial gauge transformations. 
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chirally restored deconfined phase. 4 It consists of taking two gauge fields correspond- 
ing to the two chiralities for each quark flavor in a five dimensional AdS black hole 
background. 

The action is given by two Maxwell actions for left and right gauge fields plus 
separate Chern Simons terms corresponding to separate ("consistent") anomalies 
for left and right chiral quarks. The Chern-Simons terms are, however, not unique 
but can be modified by adding total derivatives. A total derivative which enforces 
invariance under vector gauge transformations SVm = Om^v corresponds to the so- 
called Bardeen counter-term [20,21], leading to the action 

MN 1 77A P MW 



5 j^v ¥ FmA ~w/ mnFa 

+ ^mnpqr Am{f A pF A r + 3F V pF V R) ^ . (2.x) 

Since the Chern-Simons term depends explicitly on the gauge potential Am, the 
action is gauge invariant under 5 Am = Om^a only up to a boundary term. This non- 
invariance is the holographic implementation of the axial U(l) anomaly in covariant 
(Adler-Bell-Jackiw) form [29, 30] when identifying the gauge fields as holographic 
sources for the currents of global U(l) symmetries in the dual field theory. A rigorous 
string-theoretical realization of such a setup is provided for example by the Sakai- 
Sugimoto model [23,24]. As usually done in the latter, we neglect the back-reaction 
of the bulk gauge fields on the black hole geometry. 

Before we proceed, we also want to clarify our conventions concerning the e- 
tensor. We define the e tensor €mnpqr = y/—ge(MNPQR). Here we distinguish 
between the tensor and the symbol. The symbol is e(MNPQR) and normalized to 
e(r0123) = 1. 

In order to compute the field equations and the boundary action, from which we 
shall obtain the two- and three-point functions of various currents, we expand around 
fixed background gauge fields to second order in fluctuations. The gauge fields are 
written as 

A M = A M + a M , V M = V M + v M , (2.2) 

where the calligraphic letters are the background fields and the lower case letters are 
the fluctuations. 



The even simpler model considered in Ref. [26] is instead closer to a single quark flavor in a 
chirally broken phase where right and left chiralities are living on the two boundaries of a single 
branc. 
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After a little algebra, we find to first order in the fluctuations 



6S, 
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drd xy/^g < qm 



Y7 -pNM , 3K MNpQR A A y 

~2 y NJ- A -t- —e {j- np j-q R -\- j- np j- i 

.9a z 



NP) 



+VM 

+ I d 4 x 
Jd 



9v 
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(2.3) 



Qr V=9^ + 2/ee^A, J^ 



v, [ —^9K r + GK€^ X A^J X 
.9v 



From the bulk term we get the equations of motion and from the boundary terms 
we can read the expressions for the currents, 



On-shell they obey 



J» = lim —yf^T$ + 6^ upx A u T 



9v 



p\ 7 



j£ = lim —^9^ A r + 2Ke^» x A u rf x . 

r-+oo g A P 



(2.4) 

(2.5) 



V 



<W = --^ vpX (3i£F£ + F£F%) . 



(2.6) 



As expected, the vector like current is exactly conserved. Comparing with the stan- 
dard result from the one loop triangle calculation, we find k = ~^r for a dual 
strongly coupled SU(N C ) gauge theory for a massless Dirac fermion in the funda- 
mental representation. 

We emphasize that only by demanding an exact conservation law for the vector 
current can we consistently couple it to an (external) electromagnetic field. This 
leaves no ambiguity in the definitions of the above currents as the ones obtained 
by varying the action with respect to the gauge fields and which obey fl2.6| ). In 
particular, we have to keep the contributions from the Chern-Simons terms in the 
action, which are occasionally ignored in holographic calculations. 

The second order term in the expansion of the action is 



S, 



(2) 
bulk+9 



+ 



bulk 



Q-M 
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Ym 



Y7 fNM , unj MNPQR( {A -pV , fV -pA \ 

VAfJv -I- — e UNP J ~QR-t- Jnp^qr) 



(2.7) 
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where /mn is the field strength of the fluctuations. Again, the action is already in 
the form of bulk equations of motion plus boundary term. 

As gravitational background, we take the planar AdS Schwarzschild metric 

ds 2 = -f(r)dt 2 + i— + — (dx 2 + dy 2 + dz 2 ) , (2.8) 

fir) l 2 
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with / = jjj — -p-. The temperature is given in terms of the horizon by r# = L 2 7rT. 
We rescale the r coordinate such that the horizon lies at r = 1 and we also will 
set the AdS scale L — 1. Furthermore, we also rescale time and space coordinates 
accordingly. To recover the physical values of frequency and momentum we thus 
have to do replace (u, k) — > (u)/(nT), k/(7iT)). 

The background gauge fields are 

Ao(r) = $ W = «-^, (2-9) 

V (r) = *(r) = v - \ . (2.10) 

We need to relate the integration constants a, (3, 7, v to physical observables now 5 . It 
is often stated in the literature that one needs to choose a gauge in which the fields 
in Eqs. ( |2.9| ) and ( |2.10| ) vanish on the horizon in order to make AmA m and VmV m 



well defined there. This is, however, not a physical constraint. After all, the value 
of a gauge field has no intrinsic meaning. 

We will, instead, define the chemical potentials of the global U(l) symmetries 
as the potential difference between the horizon and the boundary [28]. This can be 
expressed as the integrated radial electric flux between horizon and boundary and is 
therefore a manifestly gauge invariant quantity, 

/■7\B=00 

fi = / d r A dr = MB) ~ MH) , (2.11) 

Jr H =l 

where A stands for a generic gauge potential. The variation of the chemical potential 
can be thought of as either being a variation of the gauge potential on the boundary, 
the horizon or an arbitrary combination thereof. However, by a gauge transformation 
we can always think of 5^i to result from a variation that vanishes on the horizon. 
Then 5/j, is just a special case of the general gauge field variation (|2.3|), if we interpret 
a M , v^ to be variations of the background fields. We see, therefore, that this definition 
automatically reproduces P = (Q) where Q is the integrated charge density J°. In 
general, a variation of ji is different from a variation with respect to the vector field. 
A variation in /x changes the ground state, S/i : \Q) — > \Q + SQ), whereas S/5A 
inserts the operator J° into correlation functions. 



5 Note that these integration constants with respect to the radial integration are independent of 
(t,x,y,z)\ 
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We can think of Q2-11 ) as the difference of energy in the system with a unit of 



charge at the boundary and a unit of charge at the horizon. This is the cost of energy 
to add a unit charge to the system and by definition represents thus the chemical 
potential. By the definition (j2.ll ), the integration constants (3 and 7 are thus fixed 



to 

P = /i 5 , (2.12) 

7 = fi, (2.13) 

where \x is the chemical potential of the vector symmetry and ^5 the chemical po- 
tential of the axial U(l). The constants a and v we take to be arbitrary and we will 
eventually consider them as sources for insertions of the operators J° and J® at zero 
momentum. Due to our choice of coordinates, the physical value of the chemical 
potentials is recovered by \i — > ttT/i. 

We can now compute the charges present in the system from the zero components 



of the currents ( f2.4| ) , 



,0 2 7 



J u = -i, (2.14) 

9v 

Jl = W (2 . 15 ) 

9a 

Note that this is, in fact, the standard holographic definition in the grand canonical 
ensemble. Often the gauge choice A (H) = is imposed from the outset and that 
fixes the integration constants a and v to take the values of the chemical potentials. 
It is important to realize that without a Chern-Simons term, the action for a gauge 
field in the bulk depends only on the field strengths and is, therefore, independent of 
constant boundary values of the gauge field. The action does, of course, depend on 
the physically relevant and gauge invariant difference of the potential between the 
horizon and the boundary. For our particular model, the choice of the Chern-Simons 
term results, however, also in an explicit dependence on the integration constant 
a. It is crucial to keep in mind that a is a priori unrelated to the chiral chemical 
potential but plays the role of the source for the operator J° at zero momentum. 

For the fluctuations we choose the gauge a r = 0. We take the fluctuations 
to be of plane wave form with frequency u and momentum k in x-direction. The 
relevant polarizations are then the y- and ^-components, i.e. the transverse gauge 
field fluctuations. The equations of motion are 

// ,f In/ (u 2 r 2 -fk 2 ) Ylmglk _, T , , n 

v'l + ( J j + -H + /v / J Vi + —JZ-^&Vj + *' % ) = , (2.16) 

// //' I, (u 2 r 2 -fk 2 ) YliKglk , . . N , 

a'l + { J j + -K + l p J } a % + —JL^evpOi + Vvj) = . (2.17) 
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The indices (i,j) G {y, z} and a prime denotes differentiation with respect to the 
radial coordinate r. The two-dimensional epsilon symbol is e yz = 1. 

There is also a longitudinal sector of gauge field equations. They receive no 
contribution from the Chern-Simons term and so are uninteresting for our purposes. 

The boundary action in Fourier space in the relevant transversal sector is 



S (2) 



dk 



rf 1 1 • 

— T(— O-kK)' + —V-kKY) ~ 2ikKe tJ a (a l _ k ai + 3vl fc u£) 



2^1 



<# 



(2.18) 



As anticipated, the second order boundary action depends on the boundary value of 
the axial gauge field but not on the boundary value of the vector gauge field. 

From this we can compute the holographic Green function. The way to do this 
is to compute four linearly independent solutions that satisfy in-falling boundary 
conditions on the horizon [31,32]. At the AdS boundary we require that the first 
solution asymptotes to the vector (v y ,v z ,a y ,a z ) = (1,0,0,0), the second solution to 
the vector (0,1,0,0) and so on. We can therefore build up a matrix of solutions 
Fk 7 j( r ) where each column corresponds to one of these solutions [33]. Given a set of 
boundary fields a\ (k), v\ (k), which we collectively arrange in the vector ^'^(k), 
the bulk solution corresponding to these boundary fields is 



v\Kr) = Flj^\k). 



(2.19) 



Here, F is the (matrix valued) bulk-to-boundary propagator for the system of coupled 
differential equations. 

The holographic Green function is then given by 



G /J = -2 1im(A /L (F A; L J ) / 



nj) 



(2.20) 



The matrices A and B can be read off from the boundary action as 




-2inka ( %3 

€ij 



(2.21) 



(notice that F becomes the unit matrix at the boundary). 

We are interested here only in the zero frequency limit and the first order in an 
expansion in the momentum k. In this limit, the differential equations can be solved 
explicitly. To this order the matrix bulk-to-boundary propagator is 



I 1 -9A^g(r) -9vl^g(r) \ 

9 2 A^9(r) 1 9v^9(r) 

-gyfig(r) 1 -g\^g{r) 



\ gv^gi 7 



o g 2 A ^g(r) l J 



(2.22) 



where g(r) = Giknlog(l + 1/r 2 ). We find then the holographic current two-point 
functions 

(J\p) = -I2ink(n 5 - a)ey, (2.23) 

(,r 5 Ji) = -l2iKkne tj , (2.24) 

(J l 5 ji) = -4:ink(3fi 5 - a)^. (2.25) 

Although /i, /x 5 and the boundary gauge field value a enter in very similar ways in 
this result, we need to remember their completely different physical meaning. The 
chemical potentials /i and /i 5 are gauge invariant physical state variables whereas a is 
the source for insertions of J° (0). Had we chosen the "gauge" a = [i$, we would have 
concluded (erroneously) that the two-point correlator of electric currents vanishes. 
We see now that with /i 5 introduced separately from a that this is not so. We simply 
have obtained expressions for the correlators in the physical state described by /i and 
H5 in the constant external background field a. Due to the gauge invariance of the 
action under vector gauge transformations, the constant mode v of the corresponding 
source does not appear. The physical difference between the chemical potentials and 
the gauge field values is clear now. The susceptibilities of the two-point functions 
obtained by differentiating with respect to the chemical potentials are different from 
the three-point functions obtained by differentiating with respect to the gauge field 
values. Finally, we remark that the temperature dependence drops out due to the 
opposite scaling of k and /x, /Z5. 

To compute the anomalous conductivities we therefore have to evaluate the two- 
point function for vanishing background fields v = a = 0. We obtain, in complete 
agreement with the well-known weak coupling results, 

a C ME = hm^(J^)U Q=0 = H§0 B , (2.26) 

a axial = \im^(4j% =a=0 = ^/i, (2.27) 

We are tempted to call all a's conductivities. This is, however, a slight misuse of 
language in the case of <7 55 . Formally, <7 55 measures the response due to the presence 
of an axial magnetic field B§ = V x A$. Since such fields do not exist in nature, we 
cannot measure 055 in the same way as ocme and a axial- 

Since the two-point functions in Eqs. (|2.23|) - (|2.25|) still depend on the external 



J 1 -- 


= e ct C me-B\ 


p 


CCr ax ; a ii5 , 




J5 — cr 55Bl > , 



source a, we can also obtain the three point functions in a particular kinematic 
regime. Differentiating with respect to a (and v) we find the three point functions 
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(J l (k)J j (-k)J°(0)) = 0, (2.29) 

(Jl(k)J\-k)J°(0)) = 0, (2.30) 

(Jl(k)4(-k)J°(0)) = 0, (2.31) 

(j\k)J3{-k)4m = -ik^, (2.32) 

(Jl(k)Ji(-k)J°(0)) = 0, (2.33) 

(Jl(k)4(-k)J° 5 (0)) = -ik^e^. (2.34) 

Note the independence on chemical potentials and temperature. 

Equations ( |2.32|) and ( [2.34J ) show the sensitivity of the theory to a constant 



temporal component of the axial gauge field even at zero temperature and chemical 
potentials. If the axial U(l) symmetry was exactly conserved, such a constant field 
value would be a gauge degree of freedom and the theory would be insensitive to 
it. Since this symmetry is, however, anomalous, it couples to currents through these 
three-point functions. The correlators ( |2.32| ) and ( |2.34| ) can therefore be understood 



as expressing the anomaly in the axial U(l) symmetry. 

In the next section we will check these results in vacuum at weak coupling by 
calculating the triangle diagram in the relevant kinematic regimes. 

3. Weak-coupling calculations 

An important property of the two- and three-point functions we just calculated is 
that they are independent of temperature. The three-point functions are furthermore 
independent of the chemical potentials. Therefore, the results for the three-point 
functions should coincide with correlation functions in (a chirally symmetric) vac- 
uum. At weak coupling, all the three-point functions can be obtained from a single 
1-loop Feynman integral. We only need to evaluate the diagram with two vector 
currents and one axial current. The diagram with three vector currents vanishes 
identically (due to C-parity) and the one with three axial currents can be reduced 
to the one with only one axial current by anti-commuting 75 matrices (when a reg- 
ularization is applied that permits this). Similarly, it can be seen that the diagram 
with two axial vector currents can be reduced to the one with none, which vanishes. 

When computing the three-point function, it is crucial to check the resulting 
anomalies. Gauge invariant regulators, like dimensional regularization, should yield 
the correct covariant anomaly, such that the vector currents are identically conserved. 
On the other hand, for example cutoff regularization breaks gauge invariance and 
further finite renormalizations may be needed in order to restore gauge invariance. 
In the following, we apply both dimensional and cutoff regularizations to compute 
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Figure 1: The triangle diagram. 



the three-point function and show that they give consistent results with each other 



and with Eqs. (g|2g) - qyj) . 



3.1 Triangle diagram with one axial current 

The triangle diagram, shown in Fig. [l], with one axial current and two vector currents, 
is given by 



^/J,Up 



(p,q) = (-l)(ie) 2 (ig)(i) 




i, v t+t 



(|-p)2 7 p 7 {l + qf 



r 



+(li^v,p^q). (3.1) 

The factors are a (—1) from the fermion loop, the couplings to vector and axial gauge 
fields and i for each fermion propagator. We will simply set the electric and axial 
couplings e and g to one. Evaluation of the integral with dimensional and cutoff 
regularizations is presented in some detail in Appendix [A|. 

The anomalies of the various currents coupled to the triangle diagram are ob- 
tained by contracting the three-point function above by the corresponding momenta. 
Applying dimensional regularization, we get immediately 

P,T^(p,q) = 0, (3.2) 

?,rg£(p,?) = o, ( 3 - 3 ) 



(p + q) p Kr(Pi O) 



2tt 



;PaQpe 



al3/j,u 



(3.4) 



yielding the correct Adler-Bell- Jackiw anomaly. In terms of cutoff regularization, we, 
however, find 



07T Z 

?„r&'(p, q) = ^p a q p e a ^, 

(p + i) P I co '(p»?) = -^p*m c 



r a/3fj,u 



(3.5) 
(3.6) 
(3.7) 
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In order to cancel the anomalies in the vector current, we must perform an additional 
finite renormalization by adding the Bardeen counter-term, 



.at. = c 



J d A xe^ x V^AlFj x , (3.8) 



where FK = d p V\ — d\V p . This vertex brings an additional contribution to the 
three-point function, and the full result reads 

I**' = r&'(p, q) + 2tc(p x - q x )e x ^. (3.9) 

Choosing the coefficient c of the Bardeen counter-term appropriately, c = 7^=5, we 
find the anomaly equations 

PmI^Cp, g) = 0, (3.10) 

q u T^(p,q) = 0, (3.11) 

(P + q) p r^ p (p, q) = ^P«q e a ^, (3.12) 

in full agreement with the covariant anomaly and the result from dimensional regu- 
larization. 

We next want to evaluate the triangle diagram in the special kinematic regimes 
of Eqs. (|2.29|) - (|2.34j) . Taking q = —p, corresponding to the three-point function in 



Eq. fl2.32| ), only the integrands A and B in Eqs. (|A.12|) - (|A.23|) contribute and take 
the values 1/2 and —1/2 in dimensional regularization and 1/6 and —1/6 in cutoff, 
respectively. The three-point function is then 

p-*(p p) = *ai*o p (3.13) 

2ir z 



in agreement with Eq. ( |2.32| ). Note that with cutoff regularization, | of this result 



comes from the loop diagram and | comes from the counter-term. 

Let us next take p — 0, i.e. we put zero momentum on one of the vector currents. 
The corresponding loop integral vanishes in dimensional regularization, while the 
loop contribution in cutoff regularization is precisely canceled by the contribution 
from the counter-term, 

r°^(0,-g) = 0. (3.14) 

This result is in agreement with Eq. ( |2.30D 

3.2 Triangle diagram with three axial currents 

From the same one loop integral we can also compute the correlator of three axial 
currents 6 . Since we can anti-commute the 75 and use 7! = —1, we can reduce this 



6 However, as this requires commuting the 75 with the rest of the 7 matrices, only cutoff regu- 
larization can be applied. 
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diagram to the one in Fig. (|3.1|) . The Bardeen counter-term, however, does not 
contribute this time, and we therefore find 

rf° = ^^V, (3.i5) 



just as in Eq. ( |2.34| ). The factor ~ is fixed by demanding Bose symmetry on the 
external legs. 

All other current three-point functions can be related to the triangle with three 
vector currents which is known to vanish. Therefore, we have indeed reproduced the 
holographic results in Eqs. Q2.29|) -( f2T3"4l) . 



3.3 Two-point correlator of axial currents 

To conclude the weak coupling considerations, we also compute the two-point func- 
tion of two axial currents in the background of an axial chemical potential. We simply 
can follow the analogous calculation of the chiral magnetic conductivity in [9]. The 
relevant two-point function of axial currents is 

G ^ p ) = |E/ (0^ tr b k ^ s (Qh^ s ( p + Q)] > ( 3 - 16 ) 

where S(Q) is the fermion propagator at finite temperature and density, 

S(Q) = -up — ^ ; --., (3.17) 

«7K -ip- ^575) - 11 

with to n = {2n + l)nT. Using j 5 S(Q) = —S(Q) r y5 we can square the 75 matrices to 
one and end up with the same integral as for the chiral magnetic conductivity in [9] . 
We therefore find at weak coupling that cr 55 = <tcme, coinciding with the result in 



4. Discussion and Conclusion 

We have computed two- and three-point functions of currents at finite density using 
holographic methods for a simple holographic model incorporating the axial anomaly 
of the standard model. We were able to reproduce the known weak-coupling results 
concerning the chiral magnetic effect and also found a new type of "conductivity" in 
the axial sector alone, 055. Although it can not be probed by switching on external 
fields, as a two-point function it is as well defined as ocme- It would be interesting 
to find a way of relating this anomalous conductivity to experimentally accessible 
observables. 

Previous calculations of anomalous conductivities have been able to reproduce 
the weak-coupling result for cr axia i oc fi but not oqme oc ^5 unless the contributions 



13 



from the Chern-Simons term to the chiral currents were dropped. In our calculation 
we have used the complete expressions for the currents, but of key importance was 
a clear distinction between the physical state variable, the chemical potential, and 
the external background field. The latter we viewed exclusively as a source that 
couples to an operator, whereas the chemical potential should correspond in the 
most elementary way to the cost of energy for adding a unit of charge to the system. 

It is useful to remember how a chemical potential can be introduced in field the- 
ory. One possible way is by deforming the Hamiltonian according to H —$■ H — fiQ. 
A second, usually equivalent, way is by imposing boundary conditions <p(t — i/3) = 
±exp(jU/3)0(t) on the fields along the imaginary time direction [34,35]. These meth- 
ods are equivalent as long as Q is a non-anomalous charge. Similarly, in holography 
we can introduce the chemical potential either through a boundary value of the 
temporal component of the gauge field or through the potential difference between 
boundary and horizon. Thus, for non-anomalous symmetries, the boundary value 
of the temporal gauge field can be identified with the chemical potential. Due to 
the exact gauge invariance of the action, a constant boundary value never enters in 
correlation functions. In the presence of a Chern-Simons term, however, the gauge 
symmetry is partially lost and even a constant boundary gauge field becomes observ- 



able. This can be seen explicitly from the three-point functions ( 2.32 ) and ( 2.34 ). 
Therefore, we should set the axial vector field to zero after having used it as a source 
for axial current. By defining the corresponding chemical potential as the potential 
difference between the horizon of the AdS black hole and the holographic boundary 
we are able to do so. However, the prize we have to pay is to accept singular gauge 
field configurations at the horizon. 

The fact that a gauge field that does not vanish on the horizon is not well defined 
is most easily seen in Kruskal coordinates, 

UV = - exp(47rTr*) , V/U = - exp(47rTt) , (4.1) 

where dr* = dr/f. We note that close to the horizon, r — r# rj — UV. The time 
component of a gauge field in Kruskal coordinates at the horizon is therefore 

A dt = A (r H )(^- - ^) - A' (r H )(UdV - VdU) + ■■■ (4.2) 

This is not a well-defined one- form unless Aq^h) vanishes. Although the gauge field 
is singular at the horizon, we do not believe that well defined physical observables 
are effected by this. Local gauge invariant observables, i.e. the field strengths, are 
certainly well behaved. 
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A. Evaluation of the triangle diagram 

We wish to compute the integral corresponding to the triangle diagram in Fig. [I], 

/ * / + 



■^fiup 



(p,q) 



■l)(iey(ig)(i) 




j^^Y 



+ (/i <-» v,p <H- q). 
Using Feynman parametrization the integral can be written as 

T^(p,q) = / Q/37 [tr (tbt'WVW) " tr (tsTVtVtV)] 



r 



'a/37 



-2 / dxdy 9(1 - x - y) i 



d d l N, 



a/3-y 



(2n) d (I 2 + L>) 3 ' 



where 



D = x(l — x)p 2 + 2xyp ■ q + y{l — y)q 2 



xp» - yq^ 



f 



N a/ 3j = (r - p) a rp(r + g) 7 + — 5 a(3 (r + q) 1 + 5 aj rp + 8p 7 (r - p) c 



(A.l) 

(A.2) 
(A.3) 

(A.4) 
(A.5) 

(A.6) 



Here we have already taken into account that with both dimensional and cutoff 
regularizations, the integral with odd powers of I in the numerator of the integrand 
vanishes, and the remaining tensor structure is dictated by the rotational symmetry 
of a momentum shell at fixed \l\. 

Using 

AD I A 2 



A l 5 dl 



, (P + Dr 2 
in the cutoff regular izat ion (d — 4), and 

d A ~ 2 n i 



11,8 ' ^ 
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(A.9) 
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in the dimensional regularization (d — 4 — 2e) with MS scheme, we find 



r£?(p,g) 



2tt 2 



dx cfc 9(1 - x - z) (A rcg p a + B ves q a )e 



CtflVp 



+ (CTV + Dl e ^)p a q e a ^o + (C 2 rc V + Dl e *q v )p*q^ 

with reg G {CO, DR}. The coefficients are given by 

(x — l)r 2 + yg 2 
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in the cutoff regularization, and 

(x-l)(r 2 -D) + ?/g 2 
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(1 — y)(r 2 — D) — xp 2 



A DR 
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in the dimensional regularization. 
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